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Abstract 



In this paper we study blow-up phenomena in general coupled nonlinear Schrodinger 
equations with different dispersion coefficients. We find sufficient conditions for 
blow-up and for the existence of global solutions. We discuss several applications of 
our results to heteronuclear multispecies Bose-Einstein condensates and to degen- 
erate boson-fermion mixtures. 
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ij ■ 1 Introduction 

. P.. 

The study of collapse in physical systems governed by scalar nonlinear Schro- 
dinger (NLS) equations has been a fruitful field in the last 40 years [1;2J. 
Between the physical systems for which the NLS equation is relevant to de- 
scribe concentration phenomena leading to the formation of a singularity in the 
mathematical model we can cite the optical self- focusing [51H] , Langmuir waves 
in plasmas [2] , and more recently matter waves, specifically Bose-Einstein con- 
densates [5|6|7] . 
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In other physical situations the field modelled by a Schrodinger equation is 
coupled to another field or the physical field itself is of vector nature. Examples 
of physical systems of this type in which collapse occurs are the Schrodinger- 
Poisson systems ^, the Zakharov equations [2], the generalized Zakharov 
model [9], the vector Helmholtz equation p^lllll2] . the Schrodinger-Debye 
system [13] and the Davey-Stewartson system [H]. 



However, probably the simplest and most ubiquously arising extension of the 
scalar nonlinear Schodinger equation is the vector NLS, given by the following 
equations 



idt'ijjk = --h^i'k + VkiJk+ {J2gki\^i\ Wk, (1) 

where ip^, k = 1, ...,M is a. family of M unknown complex functions defined on 
M^, D being the spatial dimensionality, A^ are the positive constant dispersion 
coefficients, Vk are external potentials acting on the different components and 
Qki are real numbers. 

Collapse in the framework of vector NLS equations has been discussed mainly 
for systems of two coupled equations (corresponding to M = 2) and Ai = 
A2 in different physical contexts such as plasma beat-wave accelerators [15] 
and incoherently coupled waves in bulk Kerr optical media [16]. The more 
general situation of an arbitrary number of components has been considered 
in Ref. [IT] for the case Afc = A in the context of copropagating nonlinear 
optical waves. The same result was rediscovered by Ghosh [18j using symmetry 
methods under similar assumptions (i.e. equal dispersions) in his study of 
the same model with applications to the dynamics of spinor Bose-Einstein 
condensates. A more abstract analysis of the same problem with only two 
components has been recently published [19j and, in the case of the absence 
of potential in [20] . 



However, several new scenarios arise in the context of collapses of matter 
waves, which do not fit into the classes of previously studied problems. The first 
one is heteronuclear multicomponent Bose-Einstein condensates [2T] where the 
atomic species have different masses contrary to the case of spinor BECs where 
different spin states of the same atomic gas are studied. In the language of Eq. 
([1]) this means that the "dispersion" coefficients for t/^i and il'2 are different, 
i.e. Ai ^ A2 [22J. 

The second field of applications of Eq. ([1]) in the theory of degenerate quantum 
gases is when a mixture of a Bose-Einstein condensate with a degenerate Fermi 
gas is studied in the mean field limit. These mixtures have received a lot of the- 
oretical attention in the last years, in particular in relation with the description 
of the observed collapse phenomena [23|24f25j . One of the models describing 
the dynamics of boson-fermion mixtures (see e.g. Refs. [26|27II28||29] ) including 
the dynamics of the fermionic component is 
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where (j)j{x,t),j 



Mp are wavefunctions describing each of the Mp 



fermions and ijj{x,t) is the mean field wavefunction for the bosonic compo- 
nent. These equations are written in adimensional units, which are related to 
the physical parameters as qbb = 4:7ih'^/mB, Qbf = ^Tih^il/niB + 1/mp), Vp 
and Vb are the magnetic trapping potentials for fermions and bosons respec- 
tively, TUB and mp are the masses of the ultracold bosonic and fermionic atoms 
respectively. The normalizations are j^n \(t)j{x,t)\'^dx = l,/ig„ \'ip{x,t)\'^dx = 
Nb- The obtention of rigorous criteria for blow-up in this systems is specially 
interesting since the large number of equations present (1 for the bosonic 
field plus Np for the fermionic fields) make the numerical analysis of realistic 
systems having thousands of fermions, unfeasible in multidimensional config- 
urations. 



Boson-Fermion mixtures as modeled by Eqs. ([2]) are also a particular case of 
our general model for multicomponent NLS systems given by Eqs. ([1]). 



Finally we would like to mention on yet another application of our model 
equations in the limit of a large number of components (formally M — >■ cxd), 
which is that of white light propagation in nonlinear media 



In this paper we study Eqs. ([T]) to obtain sufficient conditions for blow-up and 
for the existence of global in time solutions in two and three spatial dimensions. 
Our first step is to develop the theory based on the moment method (Sec. 
[2]). To do so we first define the relevant quantities for our analysis and find 
their evolution equations (Sec. 12. ip . Next, the equations related to the global 
motion of the solution are studied in Sec. 12.21 while those related to the width 
are studied in Sec. 12.31 The results on sufficient conditions for blow-up are 
then presented in Sec. [31 Finally, in Sec. H] we consider the problem of finding 
sufficient conditions for the existence of global solutions. Our conclusions are 
then summarized in Sec. O 



2 Moment equations 

2,. 1 Definition and evolution of the moments 



Our goal is to obtain a suficient condition for collapse of initial data ipk{^, t = 
0, ) satisfying equations (JT]). In what follows we will take the potentials to be 
of the form 

V,ix) = ^u;lx', (3) 

which includes the explicit forms of potentials found in realistic applications 
of NLS to the field of matter waves. 

The first part of our analysis will be based on an extended version of the 
method of moments [3T|32] for our vector system in order to obtain a virial 
type identity. The same method in a simpler form (i.e. using a smaller number 
of moments) has been used by different authors to deal with simpler situations 
(e.g. [T7|18|19ll20 ]). To apply the method we first write ipk in the modulus- 
phase representation 

^k = V^ke'''\ (4) 

and define the following functionals 



Nk= Uk dx, (5a) 

Rfc = / UkX dx, (5b) 

Pfc = / nkVipk dx = — (tpk^ipk - i^k^i'k) dx, (5c) 

Jk = J ^k {x, Vipk) dx = — I (cc, ipk'^'ipk - V'fcV^fc) dx, (5d) 

Ik= Uk x^ dx, (5e) 

Nk, = jnkn,dx, (5f) 

Kk = l {ViJk,^^k) dx, (5g) 

2n = Y.\kKk + J2collk + Y.9kiNki, (5h) 

with k, I = 1, ..., M. Hereafter integration is taken over the whole space MP . 
We will assume that all these quantities are finite for our initial data and are 
smooth as functions of time. 

After direct differentiating the previous equations with respect to time and 
some rather long but straightforward calculations we get a series of ordinary 
differential equations for the moments involving different integral quantities 



^ = 0, (6a) 

^ = A.P., (6b) 

dt 
dJi 



dt 



k 



-tUfcRfc + ^^9ki I {riiVuk - rikVni) dx, (6c) 

;/fc + XkKk - Y.9ki / nk {x, Vrii) dx, (6d) 
I •' 

— ^ = 2AfcJfc, (6e) 

— ^ = -2ujlJk -2'^gM I rik {Vipk, Vrii) dx, (6f ) 

—^ = I [XkUk {Vipk, Vrii) + XiTii {Vipi, Vuk)] dx, (6g) 
dU 

Eqs. 06ap and (16hp express respectively the conservation of the L^-norm and 
the Hamiltonian which is natural for NLS-type equations. 



2.2 Equations for the centroid of the solution 
Eqs. (I6bp and fl6c|) can be combined to get 

i d tvk 2 

Yk~dF 

where 

Ffc 



+ LuiRk = Ffc, (7) 



2^2 9ki {niVuk - UkVui) dx. (8) 

In the case of symmetric interaction matrices and homogeneous systems such 
as those studied up to now in the literature (i.e. Uk = ui, Xk = X for all 
k = 1, ..., M) these equations are a simple extension of the Ehrenfest theorem 
for multicomponent systems since one can define TZ = J2h=i -R-fc which satisfies 
the average equation TZ + uj'^TZ = 0. This equation provides a simple extension 
of the result for one component [32] and manifests the fact that the dynamics 
of the averaged center of mass behaves independently of internal forces. 

In the case of initial data with small overlapping Eqs. ([7]) can be used as 
the basis for an approximate theory for the evolution of the centroids of the 
individual wavepackets following the ideas presented in Ref. [33] . Since in this 
paper we focus our attention on collapse phenomena we will not perform such 
analysis here. 



2. 3 Equations for the width of the solution 



From Eqs. (!6el) and (]6dp we can obtain an equation for the evolution of the 
width of the individual solutions ipk which is parameter used to detect collapse 
in the framework of moment methods 



— ^ + 2ujlh = 2XkKk-2j2gkiJnk{x,Vni) dx. (9) 



Summing for all values of k and using the definitions of the momenta we find: 
Y, (^h + ^oolh) =m+{D- 2)Y,guiNu, (10) 



k ^-^^ ^ ki 

This equation is the extension of the well known virial relations for the scalar 
cubic nonlinear Schrodinger equation [2|5] for our multicomponent system 
ruled by Eqs. ([1]). 



3 SufRcient conditions for blow-up 



In this section we will use the results of the method of moments to provide 
sufficient conditions for blow-up of solutions of equations (jT]) with external 
potentials the form ^ in different scenarios. 

Before moving further we would like to remark that the Hamiltonian 7i is a 
sum of three terms H = H^-^^ + H^'^^ + H^^^ such that two of them H^^^ = 
Y, XkKkf^ and H^'^^ = Y.^l^k/'^ are always positive. The sign of the third one 

k k 

H^^^ = J29kiNki/2 depends on the properties of the matrix {gik}- Let us note 

k 

that if the matrix {gki} is such that T^*-^-* > for any Nki then collapse is not 
possible. Because of this we will be interested in the class of matrices {gki} 
which permits the negativity of the biquadratic form T-i^^' < 0. 

The question of the full description of the class of such matrices is closely 
related to the question of sign definition of the bilinear form gik^i^k but not 
equivalent to it, since in our case we have a strong cone restriction ^ > 0. We 
will not study this question in all its generality, but present one sufficient and 
easily verified characterization of a wide class of matrices for which 7i^^^ can 
become negative. 

Definition It is said that the matrix {aj^} belongs to the class 9T if there is 
a pair of indices Iq, k^ such that 

V^ioioV^koko + «iofco < 0- (11) 



— * 

One can easily check that if a matrix {gki} G ^, then there is a vector ^ 
satisfying the restriction ^j > such that Y^Qkiikii < 0. For example one can 

kl 

choose S,io = 1/Vfl'ioJo' ^ko = l/i/fi'fcofco and ^i = otherwise. 

In what follows we will assume that the matrix gik G D^ (for example, consid- 
ering the system ([2]) this means that gsF < 0) and obtain sufficient conditions 
for blow-up in that situations. In what follows we will treat the cases D = 2,3 
separately. 

3.1 Case D = 2 



Let us first study the case in which D = 2. Eq. fllOl) becomes 

Y.Y^k = ~T.^lik + m. (12) 



Let us define 



and 



J = ^4/Afc, (13) 



k 



fi2 = min (AfcCjQ . (14) 

Since A^, /fc > by definition, from ( IT2l) we get the simple evolution equation 

i + n^I = An-6'^, (15) 

where 

5' = E^(W-fi')^fc>0. (16) 

Solving flTSl) we get 

4Hl 



fc ^fc 



m 



.^(°) f^^ 



,^ , X(0) . ,^ , 47^ 
cos nt) + -^ sm(fit) + -— 



1 ^ 

- sm{n{t - t))5\t) dr. (17) 



From Eq. ([T7]) it follows that for all t < vr/fi 

J(t) < J(0) cos((]t) + ^ sin((]t) + ^ [1 - cos((^t)] . (18) 

This inequality is of the type used in the works [35|36f37j to study the stabil- 
ity and collapse phenomena in the one-component NLS equations. From this 
inequality one can conclude that if there are initial conditions such that for 



some time T* < n/Q the quantity I{T*) becomes negative, then there will 
be blow-up. The complete study of the conditions which can ensure the neg- 
ativity of the right-hand side of Eq. (1181) is a rather difficult problem. Here 
we present two (in our opinion relevant) examples of the situations when the 
analysis can be carried out. 

Example 1. Let us choose initial conditions such that 

j(0) = 0, (19) 

which can be easily done by taking the initial data to be real. Then it follows 
from Eq. (TT5]) that there exists a T* < tt/Q such that if 

n < ln^I{0), (20) 

8 

then the solution will blow-up on the time interval (0,T*). 

Example 2. If the initial data are such that 

j(0) < 0, (21) 

what always can be achieved by the appropriate choice of phases ipk of the 
initial distribution, then if 

n<\n\iio)i (22) 

then the solution blows up for times smaller or equal to T* = tt/[2Q). 

In what follows we will discuss how to find initial conditions for our equation 
([T]) with the given potentials ([3]) that can fulfill inequalities ( l20l) and f l22l) . Our 
goal is to show that the class of initial conditions satisfying fl2Ul) or fl2^ is not 
empty. The proof of this fact is done in two steps. 

Step 1. Since {gki} G Dl, then there exists a vector ^ (^j > 0) such that 
Qkiikii < 0. Choosing ip^it = 0) = \/^(t>{x) where 0(x) is an arbitrary square 
integrable function we obtain that 

^'[V'fc] < 0. (23) 

It is worth noting that such a choice of initial conditions is not a one-com,ponent 
reduction in Eqs. ([1]) because of the presence of different coefficients A^, ujj. 
and Qki (the possibility of such a reduction would imply very strong conditions 
on the coefficients A^, cj^, g^i which we do not impose). 

Step 2. Since Ti}[ipk\ can be done negative it is easy to satisfy the condition 
( 1201) . To this end we will use a rescaling argument. Let us define 



Then 



ft - \n^i{Q) = o? (n^^^ + 7^(2) - ln^i{o) - a^\n^''^ \ ] . (25) 



From the last relation it follows directly that one can always choose a suffi- 
ciently large to make (!25l) negative and thus to satisfy (!20|) . 

The case of condition (1221) is verified in a similar way. 

Remark Here we would like to make a remark concerning one very relevant 
from the physical point of view situation. Let us consider the system of equa- 
tions ([1]) in the particular case when M = 2 

idtipi = -^AiA^i + ^1^1 + r£gii H^A V^i, (26a) 



idtiJ2 = -^AsA^a + V^2^2 + (E^2« l^iH ^2- (26b) 



/=i 



From the previous considerations it follows that if the matrix {gki}k i=i belongs 
to the class 91 then the solution can blow-up in finite time. For the particular 
case of M = 2 it follows from the definition of 9T that for {gki}k i=i G 91 it is 
sufficient that 

912 < -V9u922- (27) 

This means that even in the case of repulsive intracomponent action terms 
{911, 922 > 0) blow-up phenomena can develop in finite time due to the pres- 
ence of an attractive intercomponent interaction satisfying 911922 < 9i2- More- 
over if 911,922 > and Eq. (1271) does not hold, then there are no blowing-up 
solutions of Eq. fl26l) . 



3.2 Case D = 3. 



The same type of estimates can be applied to the situation when the matrix 
{9ki} is negative definite and D = 3 since then being Nm > we can write 

i=-Y,ulik + J2 9kiNki + m, (28) 

k kl 

and thus one can obtain a similar to the two-dimensional case sufficient con- 
dition for the existence of blow-up. The same comment done in the previous 
subsection applies here, i.e. one may think of situations where a single negative 
interaction coefficient is able to drive the full system to blow-up. 



4 Global in time solutions 



In this section, assuming that the matrix {qm} has at least one negative el- 
ement, we will obtain conditions that should be imposed on the initial data 
in dimensions, D = 2, 3 to guarantee that the solution is uniformly bounded 
for all times, i.e. we will show that if the initial data are sufficiently small 
then there is no blow-up in the framework of Eqs. ([T]). The case when all the 
elements of the matrix {gki} are positive is trivial since the global solvability 
follows directly from the structure of the Hamiltonian. Moreover we will sup- 
pose that 7i > 0, since the case of 7i < implies that the initial conditions 
should be sufficiently "large" [38j (this is the consequence of the rescaling 
ipk = kipk and the structure of the Hamiltonian). 



4.1 D=3 

Let us first write the conservation of energy in the general form: 

2n = Y.Xkj\y^k?dx + Y,guj\^IJk?m^dx + Y.j'^k\uk?dx. (29) 

To obtain some useful estimations in the case of D = 3 we will now derive a 
multidimensional Gagliardo-Nirenberg inequality, namely 



EW^kMh < MC^J^Af'/'iY.W'^Ml] , (30) 

kl \ k / 

here by || • |U. we define the L^ norm : ||.|U. = (/ \v\^ dx^^ d^ stands for 
the best constant in the three dimensional scalar Gagliardo-Nirenberg inequal- 
ity [22], M is the number of components in our system and A/" = Z) llV'fclli- 
Hereafter we present a simple algebraic proof. Let us consider the Cauchy- 
Bunyakovsky inequality : 

Y. Ul^i^U's < MY. hU^- (31) 

kl '' k '' 

Now, let us apply the standard one component Gagliardo-Nirenberg inequality 
for spatial dimension D = ?> |38] 

<C^l\\^v\\%\\v\\^,, (32) 



Il4 
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where. Then from ( !3T|) we obtain that 

kl k 

Thus applying Holder inequality one more time we obtain the needed inequal- 
ity 

T.f^l^fd^^<Mcj^l^f'/^ I Y^llvMl] , (34) 

kl •' \ k / 

which corresponds to the result presented in inequality (l30l) . From this in- 
equality and the conservation of energy we can easily obtain that: 



. s3/2 

2n > Xrmn E W^Mh ' Q^MC^^IVM ^ \\^M\h , (35) 

k \ k J 

where by Xmin we denote the minimal value of all A^ and by g^, we denote the 
modulus of the minimal element of the matrix {gki}- Let us define 

U' = T.\\^i^k\\l (36) 

k 

Then 

2n > XrmnU^ " ^^AfU^ (37) 

where 7 = g^MCd. Define 

F{U) = -iVUU'^ - XminU^ + 2n>0. (38) 

The function F{U) is a cubic polynomial of U, such that F{0) = 27i > 
with two critical points Uq = and f/i = (2/3)(Amm/7wV'). Since for any 
sufficiently smooth initial the solution is locally smooth in time, then if we 
choose the initial conditions satisfying the following conditions 

• F(0) > 0, 

• FiUi) < 0, 

• U{t = 0) before the first positive root of F{U) = 0. 

then the corresponding solution is globally bounded in time. 
The first condition is satisfied automatically since 7^ > 0. As 

nUi) = -^^ + 2n, (39) 

then the second condition is equivalent to the 

Af\n\ < ^^. (40) 
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Noting that F{J2H/Xmin) > and since J2H/Xmm < Ui one can conclude 
that to satisfy the third condition it is sufficient to choose 



U{t = 0) < J^. (41) 



The conditions fj40l) and pTj) are compatible for any choice of parameters 
Afc, ujk and gik which is easily verified using the rescaling u{y) = au{(3x). 



4.2 D = 2 



Now let us analyze the case D = 2. The same ideas used to derive (|30|) allow 
to obtain that ii D = 2 then: 

Y^J^l ^f d'x < MUC^l E \\^M\l^ , (42) 



kl 



.(2) 



where CqI, is the optimal constant for Gagliardo-Nirenberg inequality in di- 
mension two. So from fl28l) we obtain that: 



271 > [Xm^n - g.MATC^l) J2 \\^M% ■ (43) 

k 

It follows from the last inequality and the constancy of the Hamiltonian, that 
if we chose initial conditions such that 

^•=Eii*.ir2.<-^. (44) 

then the solution is also globally bounded in time. 

Thus we have proved the following theorem 

Theorem 1 // the initial conditions for system of equations (Cp with the po- 
tentials of the form ^ and with the matrix {gti} having at least one negative 
element are such that 7i < and 



D=2 



^=Y.\m\%<^^^ (45) 



D=3: 

M\H\ < |-%. (46) 

27 7"^ 
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U{t = 0) < J^. (47) 



then the corresponding solution is globally bounded in time. 



5 Conclusions 



In this paper we liave studied system of coupled nonlinear Schrodinger equa- 
tions with potentials which include as particular cases models of heteronu- 
clear multicomponent Bose-Einstein condensates and mathematical models 
for Boson-Fermion mixtures. For these models we have stabilished different 
types of results on blow-up. 

First we have obtained sufficient conditions for blow-up using an extension 
of the method of moments valid for two or three spatial dimensions. In the 
most interesting situations where the coefficients of the interaction matrix 
are not all of the same sign our criterion allows to predict that there will be 
situations when the system can be driven to blowup either by a single unstable 
component or because of negative cross-interaction coefficients. 

As a second class of results we have found precise sufficient conditions for the 
initial data ensuring the existence of global solutions which makes explicit the 
intuitive idea that small enough data should be non-collapsing. 

Acknowledgements 

We want to acknowledge Vladimir Konotop (U. of Lisbon) for many fruitful 
discussions on this topic. This work has been supported by grants: FIS2006- 
04190 (Ministerio de Educacion y Ciencia, Spain), PCI08-0093 (Consejeria 
de Educacion y Ciencia de la Junta de Comunidades de Castilla-La Man- 
cha, Spain). V. Prytula is supported by grant : AP2005-4528 (Ministerio de 
Educacion y Ciencia, Spain) 



References 



[1] L. Berge, Wave collapse in Physics: Principles and applications to light and 
plasma waves, Phys. Rep. 303, 259-370 (1998). 

[2] C. Sulem and P. Sulem, "The nonlinear Schrodinger equation: Self-focusing and 
wave collapse". Springer, Berlin (2000). 

[3] V. 1. Bespalov, and V. 1. Talanov, ZhETF Pis'ma 11, 307 (1966). 



13 



[4] P.L. Kelley, Self focusing of optical beams, Phys. Rev. Lett. 15, 1005 (1965). 

[5] C. C. Bradley, C. A. Sackett, J. J.Tollett, and R. G. Hulet, Evidence of Bose- 
Einstein Condensation in an Atomic Gas with Attractive Interactions, Phys. 
Rev. Lett. 75, 1687 (1995). 

[6] C. C. Bradley, C. A. Sackett, and R. G. Hulet, Bose-Einstein Condensation of 
Lithium: Observation of Limited Condensate Number, Phys. Rev. Lett. 78, 985 
(1997). 

[7] E. A. Donley, N. R. Claussen, S. L. Cornish, J. L. Roberts, E. A. Cornell, and C. 
E. Wieman, Dynamics of collapsing and exploding Bose-Einstein condensates, 
Nature 412, 295 (2001). 

[8] F. Castella, L^-solutions to the Schrdinger-Poisson System : Existence, 
Uniqueness, Time Behaviour, and Smoothing Effects, Math. Mod. Meth. Appl. 
Sci., 8, 1051-1083 (1997). 

[9] V. Prytula, V. V. Konotop, V. M. Perez-Garcia, and V. Vekslerchik, Collapse in 
boson-fermion mixtures with all-repulsive interactions, Phys. Rev. A 76, 011803 
(2007). 

[10] G. Fibich and B. Ilan, Deterministic vectorial effects lead to multiple 
filamentation. Opt. Lett. 26, 840-842 (2001). 

[11] G. Fibich and B. Ilan, Vectorial and random effects in self-focusing and in 
multiple filamentation, Physica D 157, 112-146 (2001). 

[12] G. Fibich and B. Ilan, Multiple filamentation of circularly polarized beams, 
Phys. Rev. Lett. 89, 013901 (2002). 

[13] C. Besse and B. Bidegaray, Numerical study of self-focusing solutions to the 
Schrdinger-Debye system. Math. Model. Numer. Anal., 35, 35-55 (2001). 

[14] C. Besse and C.H. Bruneau, Numerical study of elliptic-hyperbolic Davey- 
Stewartson system: dromions simulation and blow-up. Math. Mod. and Meth. 
in Appl. Sciences, 8, 1363 (1998). 

[15] C. J. McKinstrie, and D. A. Rusell, Nonlinear focusing of Coupled waves, Phys. 
Rev. Lett. 61, 2929 (1988). 

[16] O. Bang, L. Berge, and J. J. Rasmussen, Fusion, collapse and stationary bound 
states of incoherently coupled waves in bulk media, Phys. Rev. E, 59, 4600 
(1998). 

[17] L. Berge, Coalescence and instability of copropagating waves, Phys. Rev. E, 58, 
6606 (1998). 

[18] P. K. Gosh, Exact results on the dynamics of a multicomponent Bose-Einstein 
condensate Phys. Rev. A 65, 053601 (2002). 

[19] A. Pomponio, Coupled nonlinear Schrodinger systems with potentials, J. Diff. 
Equat. 227, 258-281 (2006). 



14 



[20] L. Fanelli and E. Montefusco, On the blow-up threshold for wave weakly coupled 
nonlinear Schrodinger eqations, J. Phys. A: Math. Theor. 40, 14139 (2007). 

[21] G. Modugno, M. Modugno, F. Riboli, G. Roati, and M. Inguscio, Two Atomic 
Species Superfluid, Phys. Rev. Lett. 89, 190404 (2002). 

[22] V. M. Perez-Garcia, and J. Belmonte, Symbiotic solitons in heteronuclear 
multicomponent Bose-Einstein condensates, Phys. Rev. A 72, 033620 (2005). 

[23] R. Roth, Structure and stability of trapped atomic boson-fermion mixtures, 
Phys. Rev. A 66 013614 (2002). 

[24] T. Miyakawa, T. Suzuki, and H. Yabu, Induced instability for boson-fermion 
mixed condensates of alkali-metal atoms due to the attractive boson-fermion 
interaction, Phys. Rev. A 64, 033611 (2001). 

[25] M. Modugno, F. Ferlaino, F. Riboli, G. Roati, M. Modugno, and M. Inguscio, 
Mean-field analysis of the stability of a K-Rb Fermi-Bose mixture, Phys. Rev. 
A 68, 011601 (2003). 

[26] T. Karpiuk, M. Brewczyk, S. Ospelkaus-Schwarzer, K. Bongs, M. Gajda, and 
K. Rzewski, Soliton Trains in Bose-Fermi Mixtures, Phys. Rev. Lett. 93, 100401 
(2004). 

[27] M. Salerno, Matter wave quantum dots and anti-dots in ultracold atomic Bose- 
Fermi mixtures, Phys. Rev. A 72, 063602 (2005). 

[28] T. Maruyama, H. Yabu, T. Suzuki, Monopole oscillations and dampings in 
Boson and Fermion mixture in the time-dependent Gross-Pitaevskii and Vlasov 
equations, Phys. Rev. A 75, 13609 (2005). 

[29] J. Belmonte-Beitia, V. M. Perez-Garcia, and V. Vekslerchik, Modulational 
instability, solitons and periodic waves in a model of quantum degenerate 
bosonfermion mixtures. Chaos, Solitons and Fractals 32, 1268 (2007). 

[30] H. Buljan, M. Segev, M. Soljacic, N. K. Efremidis, and D. N. Christodoulides, 
White-light solitons. Opt. Lett. 14, 1239 (2003). 

[31] J. J. Garcia-Ripoll, V. M. Perez-Garcia, and P. Torres, Extended Parametric 
Resonances in Nonlinear Schrodinger Systems, Phys. Rev. Lett. 83, 1715 (1999). 

[32] V. M. Perez-Garcia, P. Torres, and G. D. Montesinos, The method of moments 
for Nonlinear Schrodinger Equations: Theory and Applications, SIAM J. Appl. 
Math. 67, 990-1015 (2007). 

[33] V. M. Perez-Garcia, and V. Vekslerchik, Soliton molecules in trapped vector 
nonlinear Schrodinger systems, Phys. Rev. E 67, 061804 (2003). 

[34] T.-C. Lin and J. Wei, Symbiotic bright solitary wave solutions of coupled 
nonlinear Schrodinger equations, Nonlinearity 19, 2755-2773 (2006). 

[35] L. Berge, Self-focusing dynamics of nonlinear waves in media with parabolic- 
type inhomogeneities, Phys. Plasmas, 4, 1227 (1997). 



15 



[36] L. Berge, Soliton stability versus collapse, Phys. Rev. E, 62, R3071 (2000) 

[37] L. Berge, et al., Phys. Rev. A, 62, 023601 (2000) 

[38] T. Cazenave, An Introduction to Nonlinear Schrodinger Equations, Instituto de 
Matematica, Rio de Janeiro (1996). 

[39] M, Del Pino and J Dolbeaut, Best constants for Gagliardo-Nirenberg 
inequalities and applications to nonlinear diffusions, J. Math. Pures Appl. 81, 
847-875 (2002). 



16 



